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Abstract—A one-dimensional numerical calculation has been performed on a reciprocating heat 
engine with super-adiabatic combustion in porous media. The system consists of two pistons and a 
thin porous medium in a cylinder; one being a displacer piston and the other a power piston. These 
create reciprocating motions with a phase relation, similar to those of an ordinary Stirling engine. 
By means of the reciprocating flow system, the residual combustion gas enthalpy is effectively 
regenerated to induce enthalpy increase in the mixture through the porous medium, which provides 
heat storage. The results show that due to heat recirculation, the maximum temperature is higher 
than the theoretical one in free space. Furthermore, using a porous medium with a high absorption 
coefficient, since the thermal diffusion by radiation decreases in the porous medium, the maximum 
temperature increases, resulting in the extension of the flammability limit for the heat engine. © 1997 
Elsevier Science Ltd. 
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NOMENCLATURE 

A = Frequency factor for combustion (s~') 

A, = Surface area of an equivalent particle forming porous medium (m 2 ) 

C, = Specific heat at constant volume (kJ kg -1 K _l ) 

C, — Specific heat of solid phase (kJ kg- 1 K. - 1 ) 
d = Diameter of cylinder (m) 

d, - Equivalent diameter of a particle forming porous medium (m) 

D = Diffusion coefficient (m 2 s" 1 ) 

Da = Darcy number (equation (13)) 

E = Activation energy (kJ mol 1 ) 

K, = Heat of combustion (kJ kg - ') 

He = Dimensionless heat of combustion (equation (13)) 

7/ioo = Dimensionless over-all heat transmission (equation (13)) 

HR = Dimensionless heat release rate (equation (13)) 
ke q = Permeability of porous medium (m 2 ) 

K = Overall heat transmission coefficient for heat loss (Wm -2 K.-') 

L, = Dimensionless stroke (equation 13)) 

A/ = Dimensionless heat transfer coefficient between two phases (equation (13)) 
n = Number of moles of working gas 

n, = Number density of equivalent particle forming porous medium (m -3 ) 

Nr = Radiation-conduction parameter (equation (13)) 
p = Pressure (Pa) 
r = Crank radius (m) 

r c = Thermal conductivity ratio (equation (13)) 

Re = Reynolds number (equation (13)) 

R — Gas constant (kJ mol" 1 K. - ') 
l = Time (s) 

T m Temperature (K) 
u = Velocity (ms -1 ) 

W = Reaction rate (kg m- 5 s _l ) 
x = Coordinate system (m) 

Xe = Thickness of porous medium (m) 

X = Dimensionless distance (equation (13)) 

Y = Product mole fraction 


•Corresponding author. 
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Greek symbols 

a, = Heat transfer coefficient between gas and porous medium (Wm _1 K~') 

T = dimensionless heat capacity ratio between two phases 

£1 = Emissivity of porous medium 

£2 = Emissivity of wall 

8 = Dimensionless temperature 

8 C = Crank angle (degree) 

k = Absorption coefficient (nr 1 ) 

A = Thermal conductivity of gas phase (Wm~' K' 1 ) 

A, = Thermal conductivity of porous medium (Wm _l K _l ) 

p = Viscosity (Pas~') 

p = Density of gas phase (kgm -3 ) 

p, = Density of porous medium (kg m -3 ) 

a = Stefann-Boltzmann constant (Wm -3 K -4 ) 

tp = Phase relation angle between two pistons (degree) 

cu = Angular velocity (s ') 


INTRODUCTION 

A chemical energy, such as a fossil fuel, is usually converted into mechanical and/or electric powers 
through combustion. In this case, though the combustion efficiency is close to 100%, the thermal 
efficiency and the effectiveness are not so high because of the irreversible process in combustion. 
On the other hand, it is well known that both the thermal efficiency and the effectiveness increase 
with temperature of the system. Consequently, excess enthalpy burning is a useful method for 
decreasing irreversibility of combustion and creating a high temperature flame, which will be 
successfully applied to a new energy conversion method. 

Recently, super-adiabatic combustion with reciprocating flow in a porous medium was 
investigated in an experiment [1] and a numerical calculation [2], In this system, a combustible gas 
with an extremely low heat content flows into a porous medium where the flow direction reverses 
regularly. Using the reciprocating flow system, the combustion gas enthalpy is effectively 
regenerated to induce an enthalpy increase in the combustible gas through the porous medium, 
which provides heat storage. The results revealed that the flame temperature in the porous medium 
is 13 times higher than the theoretical flame temperature of an ordinary flame in free space. 
Moreover, the heating value of the combustible gas is equivalent to a temperature increase of only 
50 K. On the basis of this novel combustion technology, a new heat engine has been proposed [3]. 
Through a numerical calculation, it is revealed that the thermal efficiency of the heat engine is 
higher than the conventional Otto and diesel cycles under the condition of the low compression 
ratio. 

The present study shows the combustion process and thermal efficiency of the reciprocating heat 
engine with the super-adiabatic combustion and shows the effects of the heat loss and radiative 
transfer to the thermal efficiency and flammability limit of the heat engine. 

CYCLE OF SUPER-ADIABATIC COMBUSTION ENGINE 

Figure 1 shows a schematic diagram of the cycle of a reciprocating heat engine capable of 
super-adiabatic combustion in a porous medium. It consists of two pistons and a thin porous 
medium in a cylinder. One is a displacer piston and the other a power piston. These create 
reciprocating motions with a phase relation, similar to those of an ordinary Stirling engine. In the 
space between the displacer piston and the porous medium, first, a fresh mixture displaces the 
combustion gas in the cylinder (i.e. scavenging stroke (Fig. 1 ©)). Then, it is compressed, similarly 
to conventional two-stroke internal combustion engines when the power piston is located in the 
vicinity of the porous medium (Fig. 1©). After the compression stroke, both pistons move in the 
same direction with almost the same velocity. As a result, the compressed mixture is preheated and 
burns in the porous medium; that is, constant-volume combustion is achieved (Fig. 1 ©). The solid 
line in the porous medium presents a schematic profile of the temperature for both the gas and 
solid phases. Combustion heat is converted into mechanical power in the expansion stroke when 
the displacer piston remains in the vicinity of the porous medium (Fig. @). On reversal, the 
residual enthalpy of combustion gas is effectively accumulated by the porous medium (Fig. 1©). 
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By means of the reciprocating motions, the residual combustion gas enthalpy is effectively 
regenerated to induce enthalpy increase in the compressed mixture through the porous medium, 
which provides heat storage. The cycle proposed here is similar to that of an internal combustion 
Stirling cycle [4-6], in which liquid fuels are directly injected into the space between the porous 
medium and the power piston and burn, similarly to conventional diesel cycles. In the present cycle, 
since constant-volume combustion is achieved, similarly to the Otto cycle, a thermal efficiency 
higher than that of the previous internal combustion Stirling cycle is expected. 

THEORETICAL ANALYSIS 

On the basis of the schematic diagram for the heat engine as shown in Fig. 1, we assumed 
one-dimensional gas flow in a cylinder after a fresh mixture displaces a combustion gas in the 
cylinder instantaneously. Both pistons are operated by means of a crank-piston system with the 
ratio of the connecting rod to the crank radius A p and with the phase relation angle i/t. The porous 
medium consists of homogeneously dispersed fine solid particles. The combustion reaction is 
described by an irreversible first-order isomerization (i.e. reactant —► product). The porous medium 
is noncatalytic. The Lewis number is unity. The physical properties are constant except for gas 
density. In the displacer piston side the surfaces of the cylinder and the piston are maintained at 
temperature T a . On the other hand, in the power piston side an overall heat transfer coefficient 
is given, including a thermal conductivity of the wall and inner and outer heat transfer coefficients; 
as a result, the surface temperature of the wall changes in a cycle. In order to clarify the effect 
of radiation in the porous medium, the diffusion approximation is applied to the analysis for 
radiative transfer. Furthermore, the radiative heat exchange between the surfaces of the porous 
medium and walls of the cylinder and piston is taken into account. The volumetric viscous force 
for the flow is expressed by Darcy’s law. 

exhaust gas 



Fig. 1. Schematic diagram of the cycle of a reciprocating heat engine capable of super-adiabatic 

combustion in a porous medium. 
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Using the above assumptions, the continuity equation and the conservation equations of 
momentum, energy and product species in the porous medium are described, as follows. 


dp d(pu) 
dt + 8x ~ u 


8u 3u dp , d 2 u u 

^ dt dx dx ^ dx 2 U 
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p = pnRT. (8) 

Herein, p and p, are the densities of gas and solid phases, respectively, p the viscosity of gas, k^ 
the permeability of the porous medium, d, the equivalent diameter of a particle, n % the number 
density of the particle, A and A s the thermal conductivities of gas and solid phases, respectively, 
h 0 the heat of combustion, a s the heat transfer coefficient around the particle, A s the surface area 
of the particle, A the frequency factor for combustion, E the activation energy, R the gas constant 
and n the number of moles of the working gas. k the absorption coefficient, a the Stefan-Boltzmann 
constant, d the diameter of the cylinder. The third terms on the right-hand sides in equations (4) 
and (6) express the heat transfer between two phases. The second term on the right-hand side in 
equation (6) expresses the radiation transferred similar to the conduction process. The fifth term 
on the right-hand side in equation (4) expresses the heat loss to the periphery. In this case, in the 
power piston side the overall heat transfer coefficient is used, while in the displacer piston side the 
inner heat transfer coefficient is used for K. The heat transfer term in equation (4), the equation 
(6) and the third term on the right-hand side in equation (2) are eliminated in calculating the 
velocity and temperature of the gas phase in the cylinder between the pistons and the porous 
medium. The boundary conditions at the high temperature side and low temperature side surfaces 
of the porous medium are given, as follows: 

—^— =---rr 4 —r 4 i x = +— (91 

’ 8x 1 I _ i ‘ ” J ’ ~ 2 K 

€l + £2 

The location of the piston surface is given by the stroke 2 r (where r is the crank radius), the ratio 
of the connecting rod to the crank radius A p and the crank angle 9 C , as follows. 

* = f + T{ 1+cos0c+ i sin ' (0c+180) }- (10) 

The boundary conditions for the gas temperature and product mole fraction are expressed as 
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T= T a , Y = 0; (11) 

for the power piston surface 

8T/8x = 0, 8Y/8x = 0. (12) 

These equations were transformed into dimensionless forms and were calculated using the 
variable finite difference method [7]. The principal dimensionless variables and parameters adopted 
here are as follows: 

X = xjx 0 , 6 = T/T 0 , U = u/rco, L, = 2r/x 0 , Re = rcox 0 p/p, H 0 = h 0 /C V T 0 , 

Da = k^/xl, HR = h 0 xl W/(C v T 0 p), M = xla s n s AJX, r c = X,/X, 

r = p s Cs/pC, Nr = kX/AgTI H ioss = 4Kd(x 0 /d) 2 /X. (13) 

For a sponge-like porous material made of cordierite, the porosity is 87.5%. Then, the values 
of dimensionless parameters are approximately T = 300, r c = 1, M = 6000 and Da = 1 x 10~ 3 , 
when x 0 = 0.01 m. The specific surface area, i.e. the product of n, and A,, of the porous medium 
is approximately 3000 m 2 /m 3 . The heat transfer coefficient a s is estimated from the heat transfer 
around a sphere, the diameter of which is equal to the mean diameter of the element forming the 
sponge-like porous medium. On the other hand, for stainless-steel piled wire meshes, the value of 
T is approximately 3000. Furthermore, the frequency factor and activation energy were estimated 
assuming that the calculated burning velocity is equal to the measured one at equivalence ratios 
of 0.53 and 1 for methane-air mixture, resulting in A = 2.6 x 10 s s~‘ and £ = 130 kJ/mol. Under 
the condition of the low velocity of the working gas considered here, the pressure drop in the porous 
medium is negligibly small for Da ^ lx 10~ 6 . The overall heat transfer coefficient K was estimated 
from the inner heat transfer coefficient 400W/m 2 K (average in a cycle [8]), the thermal 
conductivities of Si 2 N 3 ceramics (10 mm thickness) and insulator (30 mm thickness) and outer heat 
transfer coefficient, resulting in 2.34 W/'m 2 K. The emissivities of the porous medium and the wall 
are 0.2 and 0.8, respectively. The mechanical power is evaluated from the P-V diagram, i.e. the 
cylinder volume and the arithmetic mean pressure in the cylinder. 


RESULTS AND DISCUSSION 

Figures 2(a)-{f) show the profiles of the gas and solid temperatures 6, 9 S , the gas velocity U, 
the reaction rate HR and the product mole fraction Y for the crank angles 9 C of 0, 60, 120, 180, 
240 and 300 degree, respectively. Here, the dimensionless stroke L % is equal to 5, the phase relation 
angle \j/ = 60 degree, the dimensionless combustion heat H a = 3.2 and the Reynolds number 
Re = 500, which is equivalent to the speed of 300 rpm. The compression ratio is dependent of the 
stroke and the phase relation angle. In the present study, the compression ratio is 2. The 
conduction-radiation parameter Nr ( = kX/4oTI) is equal to 10. After scavenging at 9 C = 330, the 
mixture is compressed in the period over the range of crank angles from 9 C = 330-120 degree 
through the B.D.C. (bottom dead center) of the displacer piston. The mixture is preheated and 
bums in the porous medium, as shown in Fig. 2(c) and (d). The P-V diagram in Fig. 3 indicates 
that constant-volume combustion is achieved. Furthermore, as shown in Fig. 2(a), the maximum 
temperature in the porous medium is higher than the theoretical flame temperature 9 ti evaluated 
on the basis of the inlet temperature. Since the combustion continues up to the crank angle of 210 
degree, the working gas is almost isothermally expanded over the range from 9 = 180-210 degree. 
As a result, the combustion heat is converted into mechanical power, first, during the isothermal 
expansion and then during the isentropic expansion, as shown in Figs 2(d) and (e) and 3. 

With increasing conduction-radiation parameter Nr, the maximum temperature in the porous 
medium increases and the temperature profiles of gas and solid phases become steep because of 
the decrease in the effect of thermal diffusion by radiative transfer. As a result, the reaction zone 
is stabilized at the location further to the displacer piston side surface of the porous medium. 
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As shown in Fig. 3, the profiles of P-V diagram for the cases with and without heat loss are 
quite similar, though the heating values of those are not the same. The thermal efficiencies for those 
are 32% (H a = 2.7, Nr —► infinity, without heat loss) and 26% (H 0 = 3.2, Nr = 10, with heat loss), 
respectively. In the latter case, the ratio of the heat loss to the combustion heat is 38%, the most 
of which is through the cylinder wall in the displacer piston side. As a result, there is not so large 
difference between those thermal efficiencies because of a small heat loss in the power piston side; 
that is, the flowing enthalpy after expansion is divided just into the heat loss and exhausting gas 
enthalpy in the displacer piston side. 

As shown in Fig. 2, the wall temperature in the power piston side T w2 varies in a cycle. Figure 4 
shows the variation of the wall temperature in a cycle. The wall temperature has the maximum 
at the crank angle of about 180 degree, in which the reaction rate also reaches the maximum, as 
shown in Fig. 2. The maximum wall temperature in a cycle reaches the value of 5.3 higher than 
the porous medium temperature because of the compression of the working gas in the power piston 
side cylinder. However, the radiant energy emitted from the wall is absorbed by the porous medium 
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Fig. 2. Variations in profiles of gas and solid temperature, gas velocity, product mole fraction and reaction 
rate with respect to crank angle over one cycle of the super-adiabatic combustion engine. 


and regenerated to induce the enthalpy increase of mixture in the next cycle. It is considered that 
the variation of the wall temperature becomes small if the heat capacity of the wall is taken into 
account. However, the effect to the thermal efficiency is very small since there is only a small heat 
loss in the power piston side. 

Figure 5 shows the flammability limit for the super-adiabatic combustion engine, where the 
horizontal axis expresses the dimensionless heating value H 0 , while the vertical axis the 
conduction-radiation parameter Nr. With decreasing conduction-radiation parameter, the 
maximum temperature in the porous medium decreases because of the effect of the thermal 
diffusion by radiation. As a result, the flame is extinguished at the higher heating value with 
decreasing conduction-radiation parameter, as depicted by the solid line in Fig. 5. Furthermore, 
in the range between the solid and dashed lines, the maximum temperature in the porous medium 
is lower than the theoretical flame temperature in free space, though the flame is stabilized in the 
porous medium. Consequently, the super-adiabatic combustion is realized in the porous medium 
provided that the porous medium has a large absorption coefficient. 
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V 

Fig. 3. P-V diagram of super-adiabatic combustion engine. 



parameter and heating value. 


CONCLUSIONS 

A numerical calculation has been performed on a heat engine capable of super-adiabatic 
combustion in a porous medium and the following results were obtained. 

(1) On the basis of the super-adiabatic combustion with reciprocating flow in a porous medium, 
the reciprocating heat engine has a thermal efficiency of about 26% even at the flammability limit 
heating value. 

(2) The thermal efficiency is higher than that of the conventional Otto and diesel cycles for the 
compression ratio of 2. 

(3) Using a porous medium with a large conduction-radiation parameter, i.e a large absorption 
coefficient, the maximum temperature increases, resulting in the extension of the flammability limit 
for the super-adiabatic combustion engine. 
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